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Purpose ; 



To fit the log-normal tolerance distribution with quantal 
responses by maximum likelihood. 



Restrictions: 



a. The total number of objects exposed (Z n) should not exceed 61;00. 



If the number is greater, the confidence limits on the log ED, 



m, 



50' 
and the ED^q will be incorrect. Other output should be correct. 

b. At least two r. must be partial responses, i.e., for at least two 
doses 0< r- /n. < 1. 

c. Program uses normal deviates (t = Probit -5) rather than probits. 
Therefore, probit intercept, a, is coded by -5. 

d. If the magnitude of the expected t is greater than 5, the probit 
weight (w) is set equal to zero. Its actual value is less than 

0.00000771. 

e. The probit slope, b, is iterated to 3 decimal digit accuracy 
(11 binal places) or until 12 cycles of iteration are completed, 
if the slope does not converge in fewer cycles. This accuracy 
may be altered by changing the b mask in Lo + 0651+ or by extending 
the cycles by changing the maximum cycle test in Lo + 051*0. 

f . Subroutines Required: 

Routine 
Data Input No. 3 
Data Output No. 2 
Alphanumeric 
Square Root 
Log k X 
Exponential 
Antilog Printout 



Code 


Lo 


11.2 


0300 


12.1 


1300 


19.0 


1900 


15.0 


1600 


18.0 


1700 


17.0 


2300 


LED-21 


2100 



METHOD: 
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The program is stored in Lo to Lo + 0715. Data is loaded in 
Lo + 0716 to Lo + 0719 + 3N, where N is the number of data 
triplets (d,n,r). 



The method outlined in D.J. Finney: Probit Analysis , Cambridge 
University Press is used in this program. The notation used in 
the program does not always follow that of the text, but an 
appendix of notation equivalents is included at the end of this 
"writeup" . An initial estimation of a = 0, b = is used in the 
iteration. 

1 /■* -£ 
P = j J e 7 dx is approximated 

by a modification of a formula appearing in Hastings: 
Approximations for Digital Computers . The approximation is: 



p = 1 


e ^2 


vl2 Z ai u 1 
i=l 


for t < 


*J 2rr 




= 1 - [ 

Nl 


+ 2 
1 e -^- 

2tt 


5 
•42 Z a. u 1 ] 
i=l x 


for t^ 


where u = [1 + q | t 


ir 1 







a x = 0.22583681*6, a 2 = -0.252128668 
a 3 = 1.259695130, a h = -1.287822U53 
a^ = 0.9U061|6070, q = 0.23l61;l8883 
CODING INFORMATION : 

a. Storage. Program: 7 tracks plus 16 sectors (Lo to Lo + 0715) 

Intermediate and Final Results: 6320 to 632^, 6331 to 

6363. 

Data: As required - Lo + 0716 to Lo + 0719 + 3N 

b. Program Start. Lo 

c. Input. Data tape includes for each assay 

1. Compound name or number followed by a cond. stop. 

2. 0+29 (Lo + 0716) 'N'q 'q «L'-0000000» 

where N is the number of data triplets (d.,n i ,r i ) 
q^ is the q ft which doses are entered 
v 'as required) 

3 
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%. 


k 


6 


8 


10 


12 


if Zn 


0-25 


26-100 


ioi-Uoo 


U01-1600 


1601-61*00 



provided no n. exceeds 2 q n. If so, use the next higher 
even q. 

L = if not the last assay. 
= 1 for the last assay. 

3. P+qa (Lo + 0720) , d 1 'd 2 l ... 'd N « -0000000" 

k. 0+q n (Lo + 0720 + N)'n 1 , n2 l ...' n N «-0000000' 

5. 0+(q +l)(Lo + 0720 + 2N) , r 1 'r 2 ' ...' r^-OOOOOOO" 

d. Output. (See Example) 

e. Constants. Most constants are located in Lo + O6I4.J? to Lo + 0711. 

Those which may require alteration according to needs 

are: 

Maximum cycle counter (12 @ 29 ) Lo + 051*0 

Slope accuracy mask (wwwwOOOO) Lo + 0651* 

Magnitude of t 2 /2 test (12 £ @ 8) Lo + 0652* 

f. Timing. This varies according to many factors. We have found the 

average time to be 2.5 - 3»5 minutes per assay. 

g. Erasable Memory Locations. Anything outside of Lo to Lo + 0715 and 

the required subroutines. 

h. Program Stops. Lo + OO63 or Lo + 0100 indicates poor data, i.e. 

|(p-P)/Z| or |t + (p-P)/Z I > 16. 

Lo + 0200 |b I ^16. 

Lo + Q53l* - calculations completed. 

i. Format: Set CR at 1 and Tabs at 11, 21, 31 with no automatic CR 
in this range. 

j. Breakpoint 16 (Lo + 0530). 

UP - Stop after each assay, 

DOWN - No stop until the last assay is completed. 



If this is enlarged to encompass higher integers, ie. I = 13, etc, 
further values of 1 e-I must be included in Lo + 0712 ... . 



s 
m 



^ 
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Notation Changes 



D.J. Finney Program Remarks 

\ ± d. i th dose 

Y t Expected probit t = Y - 5 

a a a of program is Finney's a - $ 

Snw 2nw Sum of weights 

_ p 

Snw(x-x) [xx] Sum of weighted squared deviations 

from mean 

k N Number of data triplets (d,n,r) 

s Standard error of m (Log EDh q ) 



i l 



y_ Working probit. Coded same as a 



in program. 



Start \ 0000 




Print: 
Compound 



I 



Read from 
iape compound^ 
designation 



Data Input 3 

n, q d , q n , l, 



'Initialize* 
'djL counter, 

<Jd> d i 



5 



X i =1 °glO d i 
Replace 

<k b y x i I 




Increment 
L(d ± ) 




Set cycle 
counter b OJ 
a = 




0037 




0037 



Clear 2 nw, 
Znwx, Znwy, 

[xx],[xy],[yy] 



Initialize\ 
Kx).L(n) 
L(r) 



V 



t.=a .+b .x. 
i J J i 




Yes 



1 -«f2n 
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/Increment \ 
/L(x),L(n) \ 




n i w i & x i n 

Zn ± WjL Xi yj> 


A 


n i gi *! 
Zn ± w x y ± 


A 


n i w i y± 
& 

2n i w i 7± 


A 


2 ' 
ni Wi xi 

v & 2 
Zni Wi Xi 



4 



n ± wi Xi 
2n i w i x i 



A 



Zn t Wi 





V 




1 


Qi 


= i-Pi 


h> 


Z- 2 


U-* 


w, 

P. Q. 




011*2 



± 



X= Znwx 
Znw & 
X Znwx 



y= Znwy 

Znw 

& 
y Znwy 



A 



[xy] = 
Znwxy 
-x Znwy 



[xx] - 
Znwx 
-x Znwy 



Znwy^ 
-y Znwy 



b -&1 

j+1 [xx] 



a j + l" y - b j+l S| 
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Printout: 

#of cycles, b,a, 

x,y 

Znw, Znwx, Znwy 
[xx],[xy],[yy] 



Chi -square 
= [yy]-b[xy] 



d.f.=N-2 
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